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' Abstract. Convection plays an essential role in the emerging intensity for many stars that 

C/j , will be observed by Gaia. Convective-related surface structures affect the shape, shift, and 

asymmetry of absorption lines, the phocentric and photometric variability causing bias in Gaia 
measurements. Regarding the importance of Gaia mission and its goals, it is mandatory to 
£SJ ■ have the best models of the observed stars. 3D time-dependent hydrodynamical simulations of 

' surface convection are crucial to model the photosphere of late type stars in a very realistic way. 

I/") , These simulations are an important tool to correct the radial velocities and better estimates the 

parallaxes and photometric variability. 
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1. Introduction 

The main goal of the Gaia mission [TJ [2] is to determine high-precision astrometric parameters 
(i.e., positions, parallaxes, and proper motions) for one billion objects with apparent magnitudes 
in the range 5.6 < V < 20 and kinematic velocities of about 100 millions of stars with a precision 
of ~1 km.s" 1 up to V < 13. These data along with multi-band and multi-epoch photometric 
data will allow to reconstruct the formation history, structure, and evolution of the Galaxy. 
Convection plays a crucial role in the formation of spectral lines and deeply influences the shape, 
shift, and asymmetries of lines in late type stars which will represent most of the objects that 
will be observed by Gaia. In addition to this, granulation-related variability that is considered 
as "noise" must be quantified in order to better characterize any resulting systematic error on 
the parallax and photometric determinations. 
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Realistic modelling of stellar atmospheres is therefore crucial for a better interpretation of 
future Gaia data and, in this context, three-dimensional radiative- hydrodynamical models are 
needed for a quantitative correction of the radial velocities (few hundreds km.s -1 ) for all the 
stars observed and, in evolved stars, for the determination of the photocenter positions. 



2. Surface convection simulations and radiative transfer calculations 

We adopted here time-dependent, three-dimensional (3D), radiative- hydro dynamical (RHD) 
surface convection simulations of different stars across the Hertzsprung- Russell diagram (Tab. [I]). 
The simulations have been carried out using: 

• The box-in- a-star setup with Stagger- These models cover only a small section 
of the surface layers atop the deep convection zone, and the numerical box includes a 
number of convective cells proportional to the surface gravity. They have constant gravity, 
the lateral boundaries are periodic, and the radiation transport module relies on a Feautrier 
scheme applied to long characteristics. 

• The star-in- a-box setup with C0 5 BOLD code [3j [U [5]. These simulations cover the 
whole convective envelope of the star and have been used to model red supergiant (RSG) 
stars [3j [5] and Asymptotic giant branch stars [I] so far. The computational domain is an 
equidistant cubic grid in all directions, and the same open boundary condition is employed 
for all side of the computational box. 

The transition between the box-in-a-star and star-in-a-box occurs around logg ~ 1, when 
the influence of sphericity and the ratio between granule size versus stellar diameter become 
important; the star-in-a-box global models are then needed, but those are highly computer-time 
demanding and difficult to run so that there are only very few models available so far. 

We used the 3D pure-LTE radiative transfer code Optim3D [6] to compute spectra and 
intensity maps from the snapshots of the RHD simulations listed in Table [TJ The code takes 
into account the Doppler shifts due to convective motions. The radiative transfer equation 
is solved monochromatically using extinction coefficients pre-tabulated with the same chemical 
compositions as the RHD simulations and using the same extensive atomic and molecular opacity 
data as the latest generation of MARCS models [7J. We assumed a zero micro-turbulence since 
the velocity fields inherent in 3D models are expected to self-consistently and adequately account 
for non-thermal Doppler broadening of spectral lines. 

3. Kinematic radial velocities and correction of convective shifts 

A convenient and usual way to estimate kinematic radial velocity of a star is commonly made 
from a measurement of its spectroscopic radial velocity. Gaia is provided with a dedicated ra- 
dial velocity spectrometer (RVS, [11] ) with a resolving power at best of 11 500 centered on Ca II 
triplet in the spectral band from 8480 to 8750 A. It is well known that photospheric absorption 
lines in many stars are blueshifted as a result of convective motions in the stellar atmosphere. In 
fact, hot bright and rising (i.e., blueshifted) convective elements contribute with more photons 
than cool dark shrinking gas: as a consequence, the absorption lines appear blueshifted [12^ I13j. 
This effect is one of the systematic errors of spectroscopic measurements with RVS |14[ [15] and 
will be taken into account by the analysis pipeline of Gaia using 3D RHD simulations [16], which 
naturally account for turbulent motions and therefore do not need the use of the traditional mi- 
cro and macro turbulence used in hydrostatic models. 

We computed synthetic spectra in RVS domain using Optim3D for all the models reported 
in Tab. [1] The spectra were computed along rays of four //—angles [0.88, 0.65, 0.55, 0.34] and 
four (f)— angles [0°,90°, 180°, 270°], after which we performed a disk integration and a temporal 
average over all selected snapshots. It is important that the total time covered by the simu- 
lations is such that there are no trends in the lineshift if a subset of snapshots is used for the 
calculations. This has been verified by [17] for the simulation with T c g=4630 and log g=1.6 in 

1 www.astro.ku.dk/~aake/papers/95.ps.gz 



Table 1. 3D hydrodynamical model atmospheres used in this work. The symbols refer to Fig. 
Eld andEJ 
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Figure 1. Spectra in the Gaia-RVS domain 
for the local simulations of Tab. [TJ An offset 
has been added to the colored curves. 




o.o I . . , . i . , . . i , . . , i . . , . i , . . ^_ 

8530 B535 8540 8545 B550 8555 

A [A] 



Figure 2. Enlargements of Fig. [JJ around a 
Ca II line. 



Tab. [TJ Ramirez et al. |17] found that the sequence was long enough to ensure that the results 
were statistically significant and the comparison with observations is also very good. However, 
the other K giant models we haven't done any comparisons with observations yet. 
Figures [TJ and [2] show the results for the local simulations and Fig. [3] and [3] the results for the 
global ones. It is noticeable how the Ca II triplet (A =8498.02, 8542.09, 8662.14 A) is well visible 
in all the stars, even at low metallicity. 

Following the work of [TB], we determined the position, the depth and the width of Ca II 
triplet and of 20 Fe I lines selected to be non blended lines and to cover the spectral domain of 
the RVS [JJJ]. For this preliminary work, we fitted the absorption lines with a gaussian function. 
Figure [5] shows that Fe I are mostly blueshifted (except for few cases with low excitation 




Figure 3. Spectra in the Gaia-RVS domain Figure 4. Enlargements of Fig. [3] around a 
for the global simulations of Tab. [TJ Ca II line. 



potential). The convective lineshift (Tab[TJ AAf c i) appears stronger for high excitation lines and 
it ranges, in average, from -0.22 km.s" 1 (K giant with [Fe/H]=-3.0) to -0.75 km.s -1 (F star). 
The amplitude of the shift increases when going from K giants to the F star, as a consequence of 
the more vigorous convective motions. Moreover, the Ca II triplet lines are redshifted because 
these lines are formed in the upper part of the photosphere where the granulation pattern is 
reversed. The lineshift is particularly strong in the case of the F star (Tab[TJ AAcaii)- 
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Figure 5. Convective shifts of 20 Fe I from [19] 

as a function of excitation potentials for local Figure 6. Convective shifts of 20 Fe I from [IS] 

simulations of Tab.EJ For metalpoor stars, only as a function of excitation potentials for global 

visible lines in the spectra are reported. simulations of Tab. [TJ 



The lineshifts are not the same for all the lines in the RVS spectra. The precise amount of 
shift depends on the strength of the absorption line (and hence on the stellar metallicity), since 
different lines are formed at different atmospheric depths and thus experience different gran- 
ulation contrasts and convective velocities. In fact, Fig. [7] displays that there is a correlation 
between lineshifts and line depth indicating that weaker line (high excitation potential) are in 
general more shifted than deep lines (low excitation potential). This picture is also supported 
by Fig. In addition to this, Fig. [8] shows that there is not apparent connection with the line 
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Figure 7. Same as in Fig. Ebut as a function Figure 8. Same as in Fig. Ebut as a function 
of line depth. of i ine w idth. 



width as found by [20] for the Sun. 

RSGs show a completely different behavior (Fig. [6]) with strong redshifted Fe I up to 2.80 
km.s -1 and blueshifted Ca II up to 8 km.s -1 (Tab{T|). The reason is likely connected with the 
reversed-C shape of line bisectors observed by [21] on the prototypical RSG Betelgeuse and 
caused by a peculiar and vigorous convection but further investigations are needed and still in 
progress with RHD simulations. 

In conclusion, velocity shifts of F and K giant stars are of the order of RVS accuracy (~1 
km.s -1 ), and they are even larger for RSG stars. 

4. Photometric and photocentric variabilities and impact on parallaxes 

Massive evolved stars like RSG stars give rise to large granules comparable to the stellar radius 
in the H and K bands, and an irregular pattern in the optical region [22]. These surface inhomo- 
geneities vary with time and may affect strongly the photmetric and astrometric measurements 
of Gaia. 

Figure [9] displays the variability of a RSG during a period comparable with the Gaia mis- 
sion in the blue and red photometric bands of Gaia. The photometric system of Gaia will be 
used to characterize the star's effective temperature, surface gravity and metallicity [23J. The 
simulations show fluctuations up to 0.28 mag in the blue and up to 0.15 mag in the red bands 
over 5 years. The photometric fluctuations of RGSs are not negligible and the uncertainties on 
[Fe/H], T e g and logg for these stars with G < 15 should be revised upwards for RSGs due to 
their convective motions |24j . 

The simulated surface of RSGs (Fig. fT0|) displays high-contrast structures with spots up 
to 50 times brighter than the dark ones with strong changes over some weeks. This aspect 
is connected with the underlying granulation pattern, but also with dynamical effects such as 
shocks and waves which dominate at optical depths smaller than 1 [22]. The convective-related 
surface structures in the Gaia G band affects strongly the position of the photocenter and cause 
temporal fluctuations (Fig. left). The photocenter excursion is large, since it goes from 0.005 
to 0.3 AU over 5 years of simulation and it is on average 0.132 AU (i.e., ~ 3% of the stellar 
radius). This value corresponds to a photocenter excursion of ~10 //as for a star at 1 kpc, which 
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Figure 9. Spectral fluctuations in the blue (left) and red (right) Gaia photometric band |25] for 
a RSG simulation. The black curve is the average flux over ~5 years covered by the simulation, 
while the "grey shade" denotes the maximum and minimum fluctuations. The spectra have been 
smoothed to the Gaia photometric resolution of R ~ 50 [23]. Figures taken from 



is comparable to the Gaia accuracy of ~10/ias for stars brighter than 10. 

Chiavassa et al. [23] calculated the Gaia parallax using Gaia Object Generator v7.0, GOCQ [26] 
for RSGs with surface brightness asymmetries from RHD simulations (vj spo t) and without (w). 
Figure [TT] (right) shows that there is a systematic error of a few percent. This systematic error 
may be up to 15 times the formal error <r ro [23]. Although this error should be used to revise 
cr ro , there is little hope to be able to correct the Gaia parallaxes of RSGs from this parallax 
error, without knowing the run of the photocentric shift for each considered star. Thus, it might 
be of interest to monitor the photocentric deviations of a number of selected RSGs during the 
Gaia mission using interferometry and/or spectroscopy to obtain valuable information on the 
photocentric position and correct the resulting parallaxes using RHD predictions. 



Figure 10. Synthetic map of intensity 
(the range is [0 - 230000] erg/s/cm 2 /A) 
in the Gaia G band [25] , Figure taken 
from [23]. 




2 http://gaia-gog.cnes.fr 
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Figure 11. Ze/£: photocenter displacement as a function of time for a RSG simulation. Each 
point is a snapshot 23 days apart for a total of 5 years of simulation (comparable to the duration 
of the Gaia mission). Right: relative difference between the parallaxes computed with and 
without the photocentric motion of the left panel, as a function of the distance. Figures taken 
from [21]. 

5. Conclusions and perspectives 

3D hydrodynamical simulations of surface convection are a very important tool to account for 
convective shifts correction on the kinematic radial velocities, parallaxes and photometry. 

With the increasing computer power, it is now possible to compute 3D model grids rather 
easily on reasonable timescale. In this framework, the 3D local model grid done with Stagger- 
Code (see contribution of Remo Collet in this Volume) will be available very soon together with 
the synthetic spectra in RVS domain (before Gaia launch date, mid 2013). The grid will include 
about 130 models with T c r ranging from 4000 to 6500 K, log g from 1.5 to 5 and [Fe/H]=0.0, 
— 1.0, —2.0, —3.0. These models will not be enough to represent all the stars observed by Gaia 
and thus, following [IS], we plan to apply 3D convective corrections to spectroscopic radial 
velocities obtained with less sophisticated models and templates. 
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